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SUMMARY 


This report presents the results of an investigation concerning 
the design and testing of two-dimensional subsonic short diffusers 
employing boundary layer control using suction slots. The main goal 
was to design and test a two-dimensional Griffith diffuser using an 
improved computer design program. This program incorporates the suction 
slots in the potential flow solution. The diffuser contour was arrived 
at by prescribing the potential flow velocity along the diffuser wall 
to have a region of constant high velocity at the inlet, a region of 
constant low velocity at the exit, and a region of concentrated 
deceleration connecting the inlet and exit regions. The suction slot 
was located at the deceleration zone to control flow separation within 

this region. Side-wall suction was applied through porous side walls 
in order to maintain two-dimensional flow. Two diffusers, with area 

ratios, AR, of 3 and b, were tested. The inlet and diffuser lengths were 
1.0 in (2.5^ cm) and 6.2 in (15.7 cm) respectively. Slot suction rates 
around 8 % of the inlet flow rate were required for the diffuser with 
AR = b under metastable operating conditions, provided enough side-wall 
suction was applied. Slot suction requirements for the diffuser with 
AR = 3 were about 25% lower when a constant side-wall suction of 15% of 
the inlet flow was applied. For stable operation the required slot suction 
ranged from 28 to 30 percent for the diffuser with AR = 3 and 30 to 3b 
percent for the diffuser with AR = k. The values of the diffuser effective- 
ness n for the unseparated flow were very close to 100%. The flow emerging 
from the exit plane was uniform throughout most of the exit area. Good 
correlation was obtained between predicted and experimental values of wall 
pressure distribution and center line velocity distribution. 


1 



The secondary goal was to compare the performance of a dump (a dump 
diffuser is a sudden enlargement) and a cusp diffuser having variable area 
ratios between approximately 3 and k and having the same inlet and diffuser 
lengths as the Griffith diffuser. These diffusers were each fabricated 
with two different suction slot locations. For the dump diffuser, slots 
were provided at the inlet and exit corners; for the cusp diffuser the 
slots were located downstream of the cusp. 

For the dump diffuser (AR = 3 and k ) , the location of the suction 
slot proved to be critical regarding the diffuser effectiveness. With 
the suction slot located at the inlet corner, much higher values of 
effectiveness n were obtained than with the slot located at the exit 
corner. However, flow separation occurred in all the test runs. 

For the cusp diffuser (AR = 2.7 and 3.7), it was not possible to 
achieve unseparated flow using the available suction capability. 
Accordingly, no standing vortices were observed. The flow behaved as 
a jet emerging from the diffuser inlet section and no pressure recovery 
was obtained. 
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SECTION I 


INTRODUCTION 


1 . 1 Background 

Research on curved wall short two-dimensional diffusers employing 
boundary layer suction has been going on since 19&9 at Clemson University 
[l]'. The diffuser contour was arrived at by prescribing the potential 
flow velocity along the diffuser wall to have a region of constant high 
velocity at the inlet, a region of constant low velocity at the exit, 
and a region of concentrated deceleration connecting the inlet and exit 
regions. This contour is similar to that used on experimental high-lift 
laminar airfoils, originally suggested by A. A. Griffith [2] and is 
henceforth referred to as the Griffith diffuser. A suction slot was 
located at the deceleration zone in order to control flow separation of 
the retarded fluid. Suction through the diffuser side-walls was 
necessary in order to achieve two-dimensional flow. Test results indicated 
that a high diffuser effectiveness (98%) and a uniform exit velocity distri- 
bution were achieved. Suction rates several times higher than those pre- 
dicted theoretically were necessary to achieve these results. It was 
pointed out in [l] that the incorporation of the suction slot in the 
potential flow solution was needed. This could result in a diffuser 
requiring a lower slot suction rate than the one previously designed. 

A different approach towards designing an unseparated diffuser was 
taken by Ringleb [ 3 J • He suggested a method of preventing separation by 
means of two standing vortices located near the boundaries of the flow 

^Bracketed numbers refer to references in Section 8 
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passage. This was done by designing the contour of the flow passage in 
such a way as to allow the formation of two standing vortices with other 
flow control methods like boundary layer suction. A cusp provides the 
proper device for the standing vortices. Perkins and Hazen [A] have 
reported a successful attempt in using a short cusp diffuser. 

1 . 2 Objectives 

The primary objectives of this research effort described in this 
report were: (1) to re-examine experimentally the slot suction require- 

ment in two-dimensional Griffith diffuser with diffuser contour determined 
from the improved design program. The improved design program treats the 
suction slot geometry as an integral part of the diffuser design. A slot 
suction requirement less than ] 0 % was chosen as a limit of practical 
interest for an area ratio of I) to 1. (2) to examine experimentally the 

interaction between the side-wall suction and the slot suction requirements. 

The secondary objectives were: (1) to compare the performance, 

effectiveness and exit velocity distr i bution, of a dump diffuser (sudden 
enlargement) having the same area ratio and length to inlet width ratio as 
the Griffith diffuser at comparable suction rates. (2) to establish unsep- 
arated flow with standing vortices in a cusp diffuser having an area 
ratio comparable to the Griffith diffuser and compare its performance 
with that of the Griffith diffuser. 

1 . 3 Scope 

(1) Improve the existing design computer program which was used in the 
design of flow channels with prescribed boundary conditions (inlet, exit, 
and wall velocity distribution) to incorporate the suction slot in the 
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potential flow solution. Use the modified computer program to generate 
the contour of a Griffith diffuser. 

(2) Use the existing analysis computer program [6] to analyze the 
flow field in the diffuser geometry generated from the improved design 
program. 

(3) Use conformal mapping techniques described in Section 3-2.2 to 
arrive at the contour of a cusp diffuser having an area ratio of 4 to 1 . 

(4) Fabricate a Griffith, dump and cusp diffusers. Allowance was 
made in the design of the dump diffuser to have the suction slot location 
interchangeable between the inlet corner and the exit corner of the diffuser. 
Two suction slots were provided in the cusp diffuser, both downstream of 

the cusp. Only one suction slot was used at one time, the other was 
blocked . 

(5) Carry out experimental tests to obtain the suction requirement 
for stable and meta-stable operations, effectiveness, wall velocity 
distribution, centerline velocity distribution and exit velocity dis- 
tribution for each diffuser. Compare the measured and the analytically 
predicted velocity distributions. The range of the inlet velocity used 
was from 30 ft/sec (9 m/sec) to 260 ft/sec (79 m/sec). 
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SECTION I I 


SYMBOLS 

cross sectional area at diffuser inlet 

area ratio, exit area to inlet area 

coefficient in equation A-6 

constant in King's law 

coefficient in equation A- 7 

complex constants in equation 28 

complex potential, <j> + i \p 

used in complex numbers, defined as AT 

unit vector along x-axis 

unit vector along y-axis 

unit vector along z-axis 

outward normal unit vector 

number of segments, see equation A-6 

static pressure at any location along the diffuser wall 

stagnation pressure at end of air duct 

static pressure at inlet section 

static pressure at cusp line, see Figure 38 

static pressure at vena contracta, see Figure 38 

atmospheric pressure 

static pressure at exit section 

any point either on the body surface or in the fluid 

any point or segment on the body surface 

distance between two points p and q, see Figure A-l 

(AA. /per i meter) U. 

Reynolds number, — .j-r. —7 : : 

fluid kinematic viscosity 
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SP 


cascade spacing 


S refers to body surface; also to area on body surface; also 

length along segment in two-dimensional or length along 
stream! i ne 

s length along streamline 

U velocity 

U free stream velocity 

00 

U free stream velocity vector 

component of fluid velocity in x-direction 

U component of fluid velocity in y-direction 

U n component of fluid velocity normal to body surface 

component of fluid velocity tangential to body surface 

U. average velocity in the diffuser's inlet plane 

U e average velocity in the diffuser's exit plane 

U c average velocity at cusp line 

U t average throat velocity at vena contracta 

Uj wall velocity upstream of suction slot 

wall velocity downstream of suction slot 

U maximum velocity 

max 

Uj velocity inside the boundary layer upstream of suction slot 

V D. C . vol tage of hot wi re or hot f i 1m probe 

V q D.C. voltage of hot wire or hot film probe at U = 0 

V maxi mum D.C. vol t 
max 

W complex velocity at point p due to the segment q 

X,Y,Z cartezian coorindates 

Z complex coordinates of any point in the argand plane, z = x + iy 

Zj,Zj complex constants in equations 30 a and 30 b 
a real number in equation 27 
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Y,Y complex constants in equation 28 
6, A angles 

? complex variable defined as E, + in 

£,n real and imaginary parts of ? 

<)> potential function 

^ stream function 

p fluid dens i ty 

n diffuser effectiveness, defined by eq. (31) 
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SECTION I I I 


ANALYSIS AND DESIGN 
3.1 Two-Dimensional Griffith Diffuser 


3.1.1 Des i gn Concept 

A wall velocity distribution was prescribed as shown in Figure 1. 



Figure 1. Velocity distribution along diffuser wall. 


In region 1 and 3 an approximately constant wall velocities (with. slight 
acceleration) were prescribed with all the adverse pressure gradient con- 
centrated in a narrow region, region 2 in Figure 1. By confining the 
adverse pressure gradient to region 2, it may be possible to prevent wall 
separation in this region by applying boundary layer suction there through 
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a suction slot. This is aimed at removing the retarded fluid adjacent to 
the wall in this region. 

G. I. Taylor provided a criterion for estimating the minimum suction 
requirement across a suction slot having a concentrated adverse pressure 
gradient. Taylor's criterion as reported in [ 7 ] is given by: 


U 


1 



u 2 

(— ) 


where and are the fluid velocities outside the boundary layers just 
upstream and downstream of the slot, Uj is the velocity along the stagna- 
tion streamline just upstream of the slot, see Figure 2. 



Figure 2. Flow over a suction slot. 


In estimating the required suction rate, the boundary layer velocity 
distribution upstream of the slot is needed. Since the potential 
velocity upstream of the slot is prescribed to be nearly constant, the 
distribution may be approximated to be as that on a flat plate at con- 
stant free stream pressure. The geometry of the contoured wall can be 
obtained, assuming that the suction rate is sufficient to prevent flow 
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separation, by solving the inverse problem which is defined by prescribing 
the desired velocities at yet to be determined boundaries of the flow field. 


3.1.2 Ana lysis 

This section deals with the problem of determining the channel 
geometry for prescribed inlet, exit, and wall velocity distribution 
using two-dimensional i mcomp ress i b 1 e potential flow theory. 

Following Stanitz [8] the analysis proceeds as follows: 

Consider two-dimensional irrotational motion in the physical XY plane ) 
see Figures 3 and 4. 




Figure 3- Direction of velocity 
at a point in XY-plane. 


Figure 4. Streamlines and equi- 
potent i al 1 i nes i n 
XY-plane. 


For an element of fluid of dimensions 6s along the streamline and 6n 
normal to it, the law of conservation of mass states that: 


3 _ 

as 


(pU5n) = 0 


( 1 ) 


For an incompressible flow: 

a In U _1_ a (6n) 
3s 6n 3s 

From geometrical considerations: 

1 3(6n) 30 

— = — and 

6n 3s 3n 


1 3 (6s) o 

6s 3n 3s 


( 2 ) 


(3a) 

(3b) 


Using equations 3a and 3b into equation 2 yields: 


or 


3 In U 30_ 

3s 3n 



W 


a In U 3 £ 36 3tp 

3<j) 9s 3^ 3n 

But d \p = pUdn 
and d(j> = Uds 

Using equations 5 and 6 into equation 5, we get: 

1 3 In U 36. 

p 3$ 3\j) 


(5) 

( 6 ) 

(7) 

( 8 ) 


Equation 8 is the equation of continuity for incompressible flow in the 
- coordinate system. 

The condition of i rrotat i ona 1 i ty is: 


h (u5s) ■ 0 


(vorticity = 0) 


(9) 
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or 

1_ 96s 3 1 n U _ n 

6s 9n 9n 

Using equation 3b with equation 10 we get: 

3 In U _30_ 

3n 3s U 

or 


( 10 ) 


3 In U dip _ 3_0 3^_ 

3^ 3n 3<j> 3s 


( 11 ) 


Introducing the relationships of equations 6 and 7 into equation 11 
yields equat i on 1 2 

3 In U 39 „ (]2) 


dip 


- = 0 

a<t> 


Equation 12 is the condition of i rrotat ional i ty in the <p\p ~ coordinate 
system. 

From equations 8 and 12, equation 13 can be derived. 


3 2 In U . 3 2 In U 


(13) 


3<t> 


dip 


Thus In U satisfies the Laplace equation in the - plane and with a 
given set of boundary conditions, In U may be found over the entire 
cp-ip region. 

After equation 12 has been solved to obtain the distribution of In U 
in the transformed <j>ip - plane (for the arbitrary specified variations in 
In U with along the boundaries of constant ip), the geometry of the 
diffuser in the physical xy - plane can be determined from the resulting 
distribution of flow direction 8. The flow direction 8 can be determined 
from equation 8: 
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e 


( 1*0 


1 3 In U , 

= J p -~W~ d * 

4> 

where 41 under the integral sign indicates that the integration is taken 

along a streamline (const, ip) and where the constant of integration 
is selected to give a known value of 0 at one value of 4> along the 
stream 1 i ne. 

The distribution of the flow direction along an equ i potent i a 1 line 
(const. <p) is obtained from equation 8 which, after integration gives 




3<j> 


( 15 ) 


<t> 


where <t> under the integral sign indicates that the integration is taken 
along an equ i potent ial line ( const. <j>) and where the constant of inte- 
gration is selected to give a known value of 0 at one value of ip along 
the equ i potent ia 1 line. 

The variation in X and Y along a stream line is: 

( 16 ) 


,3X, _ ,3X 3s_x 

( d<p } , ' 3s 3<f> ; , 


Since from Figure 3 > <5X = 6s cos 0 and 6Y = 6s sin 0 and since 6<j> = U6s, 

then: 

cos 0 


3X 

9<j> 


U 


and 


Simi larly Y 


X = 


■ / 

<P 


f cos. 

u 


d(j> 


sin 8 
U 


d <}> 


( 17 ) 

( 18 ) 
( 19 ) 


The variation in X and Y along an equ ipotent i al line are given by: 


X = - 


f 


sin 0 
PU 


di^ 


(20) 



and 


V = j <H> (21) 

<}> 

The constants of integration are selected to give known values of X and 
Y at one value of <J> along a streamline or at one value of along an 
equipotential line. 

3.1.3 Incorporation of Suction Slot in the Design 

In the previous design of the Griffith diffuser which was reported in 
[l] no account was taken to the effect of slot suction on the channel 
geometry. In other words, it was assumed that suction at the walls does 
not affect the flow field or the geometry of the flow channel. This 
assumption may be justified for very small slot suction rates. The error 
thus introduced can no longer be negligible at the suction rates reported 
in [1]. An account of the way the slot suction was incorporates in the 
flow channel design is described below. 

Instead of the simple geometry used before, the wall contour near 
the slot is now as shown in Figure 5. 



=> <=> 

Figure 5. Diffuser geometry at slot location. 
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The suction flow was branched and leaving across the boundary C. As a 
consequence of this, the flow along the streamline s is brought to rest 
at a point p on the boundary just downstream of the slot. 

Boundary line C was assumed to be an equipotential line (4> - const.) 
and the flow across it was assumed uniform and parallel. From equation 15, 
it can be seen that along line C, the following conditions must be 
satisfied: 


On streamlines A and B any desired wall velocity distribution prescribed 
was as usual . 

The existence of a stagnation point p on the channel wall introduced 
a singular point; U = o and hence In U = _a> : This presented a difficult 
situation for the computer to handle. To. avoid this', Nelson [5] used 
a small portion of the potential flow solution of flow around a wedge was 
used into the finite difference mesh around the immediate vicinity of the 
wedge. In this way the stagnation point p could be isolated from the rest 
of the flow field and no infinite value encountered. 

On the <j)ij) - plane, the region around the wedge appears as in Figure 6 
below. 




The velocities at points 2, 3, and k may be expressed from the wedge flow 
solution in terms of the velocity at point 1 as follows: 


U„ = 


/M\ /T r-6s ,, 

A(j> 1 TT ' 1 


u 3 = U, 


U, = 


» - . , TT-d 

Atp /A<k .. 

’ u i 


Choosing A<}> = Aip leads to: 


U 1 " U 2 = U 3 = U 4 ” U 5 


( 22 ) 

(23) 

(2k) 

(25) 


3.1.^ Computer Design Program 

Using the analysis outlined in 3-1-2 and 3-1-3, a digital computer 
program was utilized to solve equation 13 using the Gauss~Siedel numerical 
method. The program used numerical integration to obtain the coordinates 
of the various streamlines, including the outermost streamline which 
coincide with the channel wall geometry including the suction slot. This 
assumes, of course, a negligible boundary layer thickness. The computer 
program was reported in [5] where an example was also presented to illustrate 
the details involved in preparing the input data for the design of an 
annular Griffith diffuser. A similar two-dimensional Griffith diffuser may 
be obtained merely by selecting an appropriate flag in the control card of 
the computer data deck. 

3-1-5 Computer Analysis Program • 

After the diffuser geometry for a prescribed velocity distribution was 
obtained using the design program referred to in 3-1-**, the analysis pro- 
gram was used to calculate the velocities along the diffuser wall and in 
the flow field. This provided a means of checking the diffuser contour 
generated using the design program by comparing the velocity distributions 
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calculated using the analysis program with those prescribed in the design 
program. The analysis program was developed by Hess and Smith [5] and a 
subroutine was added to the original computer program to yield the 
velocities along the diffuser wall with suction rates ranging from 1% to 
10% of the inlet flow in 1% increments. The program uses the surface - 
source - distribution method in solving channel flow problems using the 
cascade theory [6]. An outline of the theoretical analysis is given 
in Appendix A. 

3.2 Cusp Diffuser 
3.2.1 Background 

This brief exploratory study was prompted by a theory proposed by 
Ringleb (Ref. 3) on what has become known as the cusp diffuser. It is 
suggested that the snow cornice sometimes observed on the ice side of a 
mountain ridge represents a natural flow-control device which causes the 
wind on the windward side of the ridge to expand smoothly over the ridge 
onto the ice side under the influence of a vortex situated in the cavity 
behind the cusp-like edge of the cornice, see Figure 7a. This is a 
natural flow separation-control device which may prove applicable in a 
diffuser, as in Figure 7b. 



Figure 7a. Snow cornice flow. Figure 7b. Cusp diffuser. 
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The idea was used by Ringleb (Ref. 3) to build the Princeton 
University Cusp diffuser which proved, according to (Ref. k ) , to be 
workable. 

,3.2.2 Contour Generation 

The analysis makes use of conformal mapping techniques which 
transform the boundaries of the flow channel in the Z-plane into a 
simpler one in the £-plane. The flow field over the transformed 
boundaries can be readily known using a combination of sources, sinks and 
vortices. Transforming the known flow field back from the ?-plane to the 
Z-plane yields the required flow field in the channel. 

According to Riemann, for any simply connected area A in the Z-plane 
(physical plane) there exists a one-to-one conformation which maps this 
area onto the upper ?-half plane. The transformation which maps the half 
plane into the area A is then given by the function: 

Z - f(c) (26) 

where 

Z = x + i y ; and i = ^ - ! 

In transforming the?;- half plane into a divergent contour with two cusps, 
two transformations were used, these are: 

(a) l = -In [U + i ) 01 - i a ] + +f- (27) 

where a is a real number. This transforms the upper c, ~ plane into a 

divergent duct with an area ratio AR = Tr — - ■ an = 1 + — , see Figure 8. 

air a 
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equation 27 . 


For an AR = 3 
For an AR = 4 


a = 1 / 2 
a = 1/3 


(b) Z = Z* + — - + -rr-^-z (28) 

Z" - C 1 - C 

where Y, Y, C, C, are complex constants, y and C are the conjugate com- 
plex values of Y and C respectively. 

This transforms the shaded area in Figure 8b into the shaded area shown 
in Figure 9 . 



Figu re 9. 0 "|fP^| q the.sh a dgd area in Figure 8b onto the 2-plane 
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The constants in equation 28 are calculated as follows: 

(i) Select the cusp points and (Z" is the conjugate value of 

JL 

Zj) so that the cusps are slightly pointed toward the center of 
the flow passage. 

j, 

(ii) Select the constants C anc C outside the duct in the Z -plane 
so that the cusps are joining the divergent channel smoothly. 

(iii) From the condition that at the cusp points: 

■— r = 0 (29) 

dZ 


we get: 





= 0 


(30a) 


and 

1 - - ZT - 1 o — - ^ "2 = 0 (30b) 

(z, - o 2 (Zj - cr 

Equations 30a and 30b can be solved simultaneously and the constants y 
and “determined. With the constants in equation 28 now known, it is 
possible to generate the contour of the cusp diffuser. 
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SECTION IV 


APPARATUS AND INSTRUMENTATION 

4. 1 Test Faci 1 i ty 

A schematic of the test facility is shown in Figure 10. A pictorial 
view is shown in Figure II. It consists of a 10,000 CFM (4.72 m /sec) fan 
powered by a 30 HP motor, delivering air to a 20 ft. (6.1 m) long circular 
air duct having a diameter of 20 inches (50.80 cm). A shutter type 
damper at the fan inlet was used to regulate the flow rate. A 24-inch 
(61.0 cm) diameter circular plate was located at the downstream end of the 
duct. This was used as a means of connecting the diffuser to the duct. 

In order to achieve a uniform inlet air velocity to the diffusers, a 
f low-stra i ghten i ng section and a set of 4 fine screens were installed in 
the duct. The flow-straightening section was 5 ft. (1.52 m) long con- 
taining tubes of 1.25 inches (3-18 cm) diameter. The set of fine screens 
had mesh sizes of 20, 40, 50, and 100 per inch and were installed in the 
sequence of increasing mesh number in the flow direction. 

Suction was applied through the slots in the diffuser walls and 
through the side walls. This was necessary both to control separation 
and to achieve a two-dimensional flow through the diffuser. The sub- 
atmospheric suction pressure was provided by two positive displacement 
blowers each rated at 300 CFM (0.142 mVsec) at 5 inches of Hg (12.7 cm Hg) 
vacuum. The two blowers were powered by a 15 HP motor. Suction flow rate 
was controlled using a bypass line and regulating gate values. 

4.2 Gr i ff i th Di ffuser 

The computer design program [5] provided the contour of the diffuser 
(AR = 4). Figure 12 shows the coordinate and the contour resulting from 
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the design program. A computer plot of the diffuser contour is 
shown in Figpre 13. Figure 14 shows the wall velocity distribution 
prescribed to the design program with 5 % of the through flow as the 
design suction rate for the suction slots, and for AR = 4. 

The diffuser was fabricated from aluminum stock to tolerances within 
one-thousandth of an inch (.0025 cm). The side walls were fabricated 
from a sintered porous stainless steel plate. A pictorial view of the 
diffuser, attached to a circular plate, is shown in Figure 15. The 
circular plate provided a means of connecting the diffuser to the air 
duct. 

4. 3 Dump D i f fuser 

The dump diffusers used were simply sudden enlargements with area 
ratios of 3 and 4. Suction slots were provided either at the inlet or 
exit corner; Figure 16 shows two revisions of a dump diffuser with the 
slots located at either the exit corners or at the inlet corners. 

4. 4 Cusp Di f fuser 

Using the method outlined in 3.2.2, the contour was generated using 
the transformation equation 28. The constants C and y used are: 
c = 0.150 + 1 . 050 i , y = -0.134 + 0.208i 
Figure 17 shows the resulting coordinates and contour. Two suctions slots 
were provided on each wall on the downstream side fo the cusp. This 
would help to stability two trapped vortices [3]. Side-wall suction 
was also provided. The diffuser had the same inlet and diffuser length 
as both the Griffith and dump diffusers, 1 inch (2.54 cm) for inlet 
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length and 6.2 inch (15-7 cm) for diffuser length). Using these 
dimensions it was not possible to design a cusp diffuser with an area 
ratio of 4; the resulting diffuser had an area ratio of 3.7. This 
approximated the desired area ratio of 4 and may be set as 2.7 to 
approximate 3 by moving the two halves of the diffuser walls relative to 
each other in a direction normal to the air flow. Figure 18 shows a 
pictorial view of the two contoured walls of the cusp diffuser. 

4.5 Instrumentation and Measurements 

The stagnation pressure at the end of the air duct and the static 
pressure at the diffuser inlet were measured using a Mer i am- type-M i cro 
micromanometer with a resolution of 0.001 inch (0.0025 cm) of water. 

Using the Bernoulli equation and the difference of the above measured 
pressure readings the inlet velocity was determined. 

The wall pressure distribution was measured using a differential 

2 

pressure transducer with a range 0 - 1.0 psid (0 - 0.69 N/cm ) having an 
error of 0.5% of the full scale, along with a rapid scanning mechanism 
with a capacity of 48 channels. The number of wall pressure taps on each 
wall was as follows: 


Gr i ff i th Di f fuser 

17 

Dump 

Di f fuser 

18 

Cusp 

Diffuser 

17 


The output voltage from the pressure transducer was measured using a 
Digitec 269 Multimeter digital voltmeter and automatically recorded on 
paper tape for subsequent analysis. 

The exit and center line velocity distribution were measured using 
hot wire and hot film probes in conjunction with a constant temperature 
anemometer. The use of either probe depended mainly on availability. 
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The D.C. voltage output of the probe was displayed by a multi-range 
Disa-type 55 D30 digital voltmeter. The probe was positioned using a 
traversing mechanism which had a resolution of 0.001 inch (0.0025 cm) 
in each of the three perpendicular directions. 

The exit velocity distribution was obtained by making side-wall 
to side-wall velocity traverses along seven rows in the exit plane, 
see Figure 19. Along each row 19 measurements were taken which makes 
the total number of exit velocity measurements to be 133- Exit velocity 
distributions for the area ratio of 3 and 4 were obtained. In addition 
to the exit velocity distributions, axial velocity measurements were 
taken at several positions along the diffuser axis in the flow direction. 
The main through flow at the diffuser inlet section was estimated from 
the measured pressure difference between the inlet section and the plenum 
chamber at the end of the air duct. The exit velocity distribution and 
the flow through the suction passages measured by using venturi meters 
provided a check on the value of the inlet flow estimated by using the 
inlet area and the inlet velocity measurement mentioned earlier. 


25 



Schematic of test facility 









Figure 11. Pictorial view of test facility looking downstream 



Figure 12a. Contour of Griffith! Figure 12b. Coordinates of Griffith 
diffuser diffuser 


Figure 12. Coordinates and contour of Griffith diffuser 
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Figure 13. Computer plot for Griffith diffuser contour showing streamlines and equ i potenti al 
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Slots at inlet corners Slots at exit corners 


Figure 16. Sketch of dump diffuser showing two versions 
of suction slot locations 
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Figure 17a. Contour of cusp diffuser 


Figure 17b. Coordinates of cusp diffuser 
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SECTION V 


TEST CONDITIONS, PROCEDURE AND DATA REDUCTION 


5.1 Test Cond i t i ons 

To achieve the outlined project objectives, a comprehensive 
experimental investigation was carried out for the Griffith diffuser. To a 
lesser degree, experiments were conducted on both the dump and the cusp 
d i f fusers . 

Table 1 gives a summary of the range of test parameters used in the 
tests for all three types of diffusers. 


Table 1: Range of Test Parameters 


Test Parameter 

Range 

From 

To 

Inlet air velocity, ft/sec 

30 

260 

Inlet air velocity, m/sec 

9 

79 

Inlet volumetric flow rate, CFM 

0 

2500 

Inlet volumetric flot rate, m3/min 

0 5 

70 6 

Inlet Reynolds number 

6 X 10 p 

7 X 10 6 

Suction volumetric flow rate, CFM 

0 

600 

Suction volumetric flow rate, 



m^/mi n 

o,. 

17.,. 

Area ratios used 

3 a 

id 4 

Slot suction (percent of inlet 



flow rate) 

0 

50 

Side-wall suction (percent of 



inlet flow rate) 

0 

35 

Slot width, inch 

1/16 a 

nd 1 /4 

Slot wi dth , cm 

0.16 and 0.64 

Slot length, in 

8 

Slot length, cm 

20.3 

1 


2.7 and 3-7 for cusp diffuser 


For all diffuser tests, room air was admitted to the fan through the 
shutters and delivered through the air duct, straightening tubes and 
screens to the diffuser and discharged back into the room. 
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5.2 Test Procedure and Data Reduction 


Each test was conducted to obtain the following informations in 
the sequence indicated: 

(a) Find the minimum slot suction requirements for both meta- 
stable and stable operating conditions for different inlet air velocities 
and for a constant side wall suction rate. 

(b) Find the value of the diffuser effectiveness n for different 
values of inlet air velocity for the case of attached flow. 

(c) Find the diffuser wall velocity distribution by taking wall 
static pressure readings using the pressure transducer. 

(d) Find the exit velocity distribution by making horizontal 
traverses at the exit plane by taking exit velocity measurements using 
the hot wire or hot film probe. 

(e) Find diffuser axial center line velocity distribution by 
making an axial traverse. 

The following procedure and relationships for data reduction were 
used. In order to obtain the minimum slot suction requirement for 
meta-stable operating condition for a particular inlet air velocity the 
suction rate was first increased to achieve attached flow. The suction 
rate was then gradually reduced until the flow pattern in the diffuser 
changed from the attached to the separated condition as indicated by tufts 
on the exit walls. At this point the suction rate was slightly increased 
to restore the flow to the attached condition and all the necessary 
measurements were taken at the diffuser exit. Any interference 
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with the boundary layer at the slot or upstream of the slot would cause 
the flow to be separated. Therefore no probe was inserted in this region. 

The minimum slot suction requirement for the stable operating condition 
for a particular inlet air velocity was obtained by increasing the slot 
suction beyond that required for minimum meta-stable operating condition. 
With the attached flow condition now achieved, a test for stability was 
performed. This was done by physically disturbing the boundary layer at, 
and upstream, of the slot; for example by touching these areas by hand. 

This would cause an immediate boundary layer separation. If slot suction 
was enough to maintain stable operation, the flow would be restored to 
its previous attached condition once the disturbance was removed. 

Otherwise, the suction had to be increased and the test repeated until 
the condition of stability obtained. 

The diffuser effectiveness was estimated from three independent 
pressure measurements made using the micromanometer: (l) the static 

pressure at the diffuser inlet, P.; (2) the static pressure at the end 

of the air duct, P q ; (3) the static pressure at the diffuser exit, 

P . From these three pressure measurements and from the previously 
measured total suction requirement, the diffuser effectiveness n was 
calculated from: 



The diffuser wall velocity distribution was estimated from the wall 
static pressure measurements provided by the differential pressure 
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transducer. The transducer provided values for the pressure difference 
between the stagnation pressure P q at the end of the air duct and the 
pressure P at any particular location along the diffuser wal 1 . From 
these measurements the ratio of the wall velocity U at any particular 
location to the wall velocity U. at the inlet section was calculated by: 



P 

PT 


' V o i 

Both the exit velocity distribution and the center line velocity 
distribution were estimated from either constant temperature hot wire 
or hot film probe measurements using King's law: 


V 2 = + Afu, 


V i s the D.C. vol tage 


at any velocity U and is the D.C. voltage at U = 0 and A is a constant. 
Several probes were calibrated against pitot-static tubes and the results 
justified the use of King's law. The ratio of the velocity at any location 


U to the maximum velocity U is given by: 

max 3 7 

2 2 

V - V 

o 


max 


2 2 

W - V 

max o 


where V is the D.C. voltage corresponds to U 

max 3 r max 
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SECTION VI 


DISCUSSION OF RESULTS 
6.1 Griffith Diffuser 

Figure 20 shows the minimum slot suction requirement for con- 
stant values of side-wall suction rate expressed in percent of volu- 
metric inlet flow rate for the diffuser with area ratio, AR = A, and 
slot width = 1/8 inch (0.32 cm). The value of minimum slot suction 
requirement for the metastable condition and for a particular inlet 
velocity depends on the percent side-wall suction applied. At an 
inlet velocity of 100 ft/sec (30.5 m/sec) the slot suction require- 
ment drops by 50% by increasing the side-wall suction from 10% to 
25%. The amount of reduction in slot suction decreases as the 
side-wall suction continues to increase. This becomes obvious if 
one examines the effect of increasing the side-wall suction 



Figure 20. Percent-Slot suction versus inlet velocity for Griffith 
diffuser; AR = 4; slot width = 1/8 in (0.32cm) 
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from 10% to 25% in steps of 5% for an inlet velocity of 100 ft/sec 
(30.5 m/sec). An increase of side-wall suction of from 10% to 15% 
causes a reduction in slot suction of 28%, while an increase of from 
20% to 25% causes a reduction of only 5%. Due to limitations in 
suction capability, it was not possible to achieve stable operating 
condition for more than few cases. Figure 21 shows the minimum meta- 
stable slot suction requirements for the same area ratio, AR = 4, 
but with a slot width double that of Figure 20; 1 /b inch (0.64 cm). 

The same trend is observed, i.e. an increase in the percent side-wall 
suction reduces the minimum slot suction requirement to achieve meta- 
stable operating condition. In this case; slot width = 1/4 inch 
(0.64 cm); the values of the slot suction requirements seem to be 
higher than the values for the case where Slot width = 1/8 inch (0.35 cm). 



INLET VELOCITY (t/,ec 


Figure 21. Percent slot suction versus inlet velocity for Griffith 
diffuser; AR = 4; slot width. = 1/4 in (0.64 cm) 
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For a side-wall suction of 25%, the slot suction requirement is ahout 
IA% for the wide slot while for the narrow slot it is only about 8%. 

The effect of area ratio on the minimum suction requirement is shown 
in Figure 22 for two area ratios, AR = 3, and for a 15% side-wall 
suction. As would be expected, the diffuser with the smaller, area ratio; 

AR = 3; requires a smaller percent slot suction to achieve meta-stability 
than would the larger one with AR = k. The reduction in meta-stable slot 
suction requirement for the small area ratio, AR = 3, is about 25%, for 
the same side-wall suction of 15%. Using Taylor's criterion, and the 1 /7th 

power law to approximate the velocity distribution in the boundary layer, 

the percent reduction in slot suction for AR = 3 was estimated to be 18%. 

For stable flow, the required suction rate ranged from 28 to 30 percent 

for the diffuser with AR = 3. 

It is apparent that slot suction percentages within a practical 
range can be achieved at the expense of higher side-wall suction per- 
centages. It is believed that high side-wall suction was necessary 
because the boundary layer thickness at the corner between the diffuser 
walls and the side walls was several times thicker than elsewhere. 
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Referring to Figure 20.13 page 576 of Sch 1 i cht i ng 1 s Boundary Layer Theory 
[9] one notices that the boundary layer thickness at the corner of a 
noncircular duct is several times thicker than along the flat surfaces. 
Accordingly, it was more likely that separation would start at the 
junction of the diffuser walls and the side wallsi Although no 
attempt was made at applying more suction along all the corners, several 
trials were made aiming at reducing the side-wall suction requirement 
by drilling several small holes in the sintered porous side plates in the 
region where the slots meet the side plates. This would apply more 
suction to the junction between the slot and the side walls and reduce the 
chance of separation in this region. No appreciable reduction in required 
side-wall suction was achieved. 

The diffuser effectiveness, n, of 14 test runs are shown in Figure 23 
over an inlet air velocity range of 50 ft/sec (15.2 m/sec) to 180 ft/sec 
(5^-9 m/sec) for both area ratios of 3 and 4 and for slot widths of 



INLET VELOCIJT *»/••« 


Figure 23. Effectiveness q versus inlet velocity for Griffith 
diffuser; AR = k. 
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DIFFUSER EFFECTIVENESS 


1/8 inch (0.32 cm) and 1/4 inch (0.64 cm). A 100% pressure recovery 
is achieved using the suction percentages shown in Figures 20 and 21. 
Similar values of q were reported in [1]. In Figure 24 the diffuser 
effectiveness is shown for varying values of slot suction with constant 
side-wall suction for stable and meta-stable operations. At values of 
slot suction higher than the minimum required in each case, 
the effectiveness has a value of 100%. As the value of slot suction is 
reduced below the minimum, the effectiveness drops sharply from 100% to 
about 15%. This is the value of the effectiveness with separated flow. 

The wall velocity distributions as calculated from the diffuser 
wall static pressure measurements are shown in Figures 25 and 26 for 



0 4 8 12 16 20 24 28 32 

PERCENT SLOT SUCTION 


Figure 24. Effectiveness n versus percent slot suction for Griffith 
di ffuser ; AR = 4. 
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LOCAL MALL VELOCITY / INLET VELOCITY 



Figure 25- Wall velocity distribution for Griffith 
diffuser; AR = 3* 



Figure 26. 


Wall velocity distribution for Griffith 
diffuser; AR = 3. 




area ratios of 3 and 4 respectively. Also shown on the same figures 
are the wall velocity distributions prescribed to the computer design 
program with 5% slot suction. This would help to compare the trend 
rather than the values. The agreement between the measured and pre- 
scribed values of the wall velocities downstream of the suction slot 
is misleading. This is because of the difference in total suction be- 
tween the design value and that used in the test run. For example, 
referring to Figure 25, the total suction used, which is the suction 
used for the slots and the side walls, is 24% of the inlet flow rate 
while the design value is only 5%. Accordingly one would expect the 
experimental results of the wall velocities downstream of the slot 
to have values lower than those prescribed to the design program. Re- 
ferring again to Figure 25, the ratio of the exit velocity to inlet 
velocity calculated from the wall static pressure measurements is found 
to be 0.23- A simple flow balance yields a value of 0.25 for the same 
quantity which is 8% higher than the previous value (0.23). This can 
be attributed to the fact that the measured values of wall velocities 
were normalized by the inlet wall velocity which is higher than the 
average inlet velocity used in the flow balance calculations. A simi- 
lar check is made in the case of Figure 26. 

The results of the center line velocity traverses are shown in 
Figures 27 and 28 for the area ratios of 3 and 4 respectively. Also 
shown are the distributions resulting from the computer analysis pro- 
gram using values for slot suction equal to the total suction used in 
the test runs; i.e. 27% for AR = 3 (Figure 27) and 42% for AR = 4 
(Figure 28 ). The agreement between the analysis and test results is 
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Figure 27. Center line velocity distribution for 
Griffith diffuser; AR «= 3- 
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Figure 28 


Center line velocity distribution for 
Griffith, diffuser; AR = k. 





fairly good. Referring to the same figures, flow balance calculations 
yield the values of 0.24 and 0.14 for the ratio of the average exit to 
inlet velocity and for the area ratios of 3 and 4 respectively. The 
measured ratio for the center line velocities are slightly higher than 
0.24 and 0.14 as can be seen from the same figures. This can be attri- 
buted to the fact that the center line velocity at the inlet section is 
lower than the average velocity at the inlet section; the velocity dis- 
tribution exhibits a "dip" in the center portion. 

Figures 29 and 30 show exit plane velocity maps for the diffuser 
with area ratios 3 and 4 when it is operating in the meta-stable con- 
dition. The inlet air velocity was kept constant at 170 ft/sec (51-8 
m/sec) and the suction was adjusted to obtain a meta-stable operating 
condition. The contours of constant exit plant velocity exhibit a flat 
profile for the major portion of the exit plane. For both area ratios, 
a slight "dip" in exit plane velocity was observed in the area around 
the center while velocities near the diffuser and side walls are general- 
ly higher. Figure 31 shows the exit plane velocity distributions along 
the central horizontal line, C-L. The exit plane velocity distribution in 
the vertical direction obtained from the computer analysis program and the 
computer design program are shown in Figure 32, together with the measured 
values. The discrepancy between the theoretical and experimental results 
could be due to a certain degree of nonuniformity in the flow through the 
air duct which was transmitted to the diffuser inlet. 

6.2 Dump Diffuser 

It was not possible to operate the dump diffuser unseparated either 
in the meta-stable or stable condition using the suction capacity avail- 
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Figure 29. Exit plane velocity map for Griffith 

diffuser; AR - 3 (lines of constant U/U ) 
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Figure 30. Exit plane velocity map for Griffith 

diffuser; AR = k (lines of constant ll/ll ) 
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Figure 31. Exit plane velocity distribution 
along horizontal center line for 
Griffith diffuser; AR = 3 and 4. 



Figure 32. Exit plane velocity distribution 
along vertical center line for 
Griffith diffuser; AR = 3 and 4. 
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able. All measurements were taken in the separated condition. The 
area ratio was varied from 3 to if by varying the spacing between the 
two diffuser walls. Two suction slot locations were tested, one at 
the inlet corner and the other at the exit corner. The suction slots 
were 8 inches (20.3 cm) long and 1/16 inch (0 . 1 6 cm) wide. 

Figures 33 and 34 show the values of the effectiveness obtained 
for the dump diffuser with an area ratio of 3 and 4 with the suction 
slot located at either the inlet corner or the exit corner. The slot 
width used for both locations was 1/16 in (0 . 1 6 cm). The slot suction 
and side-wall suction varied from 9% to 47% and 5% to 30% respectively. 
The particular values of slot and side-wall suctions for each test 
run are shown in Table 2. Comparing the values of the effectiveness ob- 
tained with the suction slot located at the inlet and exit corners 
shows that substantially higher values were obtained when the slot was 
located at the inlet corner. 

Figures 35 and 36 show the wall velocity distribution for area 
ratios of 3 and b respectively and with the suction slots located at 
the inlet corners. The test conditions for both cases are: inlet 

velocity = 148 ft/sec (45-1 m/sec), slot suction = 33% and side-wall 
suction = 20%. 

The exit velocity distributions are shown in Figures 37 and 38 
for the area ratio of 3 and 4 respectively and for the same test condi- 
tions of Figures 35 and 36. These show a large degree of non-uniformi- 
ties in comparison with the exit velocity distribution of the Griffith 
diffuser (see Figure 30- 
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DIFFUSER EFFECTIVENESS 


INLET VELOCITY m / sec 



INLET VELOCITY ft/sec 

Figure 33- Effectiveness n versus inlet velocity for 
dump diffuser, AR = 3- 


INLET VELOCITY m/ sec 



INLET VELOCITY ft/ sec 

figure 34. Effectiveness n versus inlet velocity for 
dump diffuser, AR = 4. 
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Table 2. Test Conditions and Effectiveness of Dump Diffuser 

AREA RATIO, AR = 3 

Slot At Inlet Corner Slot At Exit Corner 


Inlet 

ft/sec 

Veloci ty 
m/sec 

Suet i on 
Slot 

% 

Side 

Wall 

Ef fect- 
i veness 
n 

Inlet 

ft/sec 

Veloci ty 
m/sec 

Suction % 
Slot Side 

Wall 

Effect- 
iveness n 

44.3 

13.5 

17 

24 

0.88 

105-0 

32.0 

34 

30 

0.23 

86.8 

26.4 

1 1 

13 

0.83 

157.0 

47.8 

22 

21 

0.11 

87.0 

26.5 

14 

15 

0.86 

203.0 

61 .9 

17 

16 

0.08 

1 16.0 

35.4 

31 

28 

0.83 

230.0 

70. 1 

15 

14 

0.08 

131-0 

39.9 

14 

9 

0.91 






131-0 

39.9 

19 

12 

0.89 






154.0 

46.9 

23 

21 

0.71 






173-0 

52.7 

14 

1 1 

0.85 






173.7 

52.9 

1 1 

9 

0.82 






205.0 

62.5 

17 

16 

0.61 






232.0 

70.7 

9 

4 

0.75 






232.0 

70.7 

11 

5 

0.75 









AREA RATIO, AR = 4 





Slot At Inlet Corner 



Slot 

At Exi t 

Corner 












88.0 

26.8 

17 

17 

0.85 

104.0 

31.7 

47 

29 

0.12 

132.0 

40.2 

32 

14 

0.89 

210.0 

64.0 

23 

14 

0. 16 

148.0 

45.1 

33 

20 

0.77 

230.0 

70. 1 

21 

13 

0.11 

175.0 

53.3 

22 

10 

0.84 

250.0 

76.2 

19 

12 

0.08 

178.0 

54.2 

16 

1 1 

0.80 






222.0 

67.7 

22 

14 

0.65 






257.0 

78.3 

19 

12 

0.57 
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dump diffuser; AR = / f. 
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Figure 36. Wall velocity distribution for 
dump diffuser; AR = k. 


50 





EXIT VELOCITY RATIO, LOCAL/ MAXIMUM EXIT VELOCITY RATIO, LOCAL/ MAXIMUM 



AREA RATIO * 3 
SLOT WIDTH * 1 IN 3 -16 cm 
16 

SLOT SUCTION ■= 33 % 

SIDE WALL SUCTION = 20 % 

INLET VELOCITY * 148 0 ft/sec 
= 451 m / 86C 

O CENTER LINE 
A I 5 IN (3-8 cm) FROM TOP 
* I 5 IN (3 8 cm) FROM BOTTOM 



HORIZONTAL DISTANCE (x) / WIDTH(X) 


Figure 37. Exit plane velocity distribution along horizontal lines 
for dump diffuser; AR = 3* 



AREA RATIO = 4 

SLOT WIDTH = 1 IN » -16 cm 

SLOT SUCTION = 33 % 

SIDE WALL SUCTION 3 20% 

INLET VELOCITY » 148 0 ft /sec 
° 45 I m / sec 

O CENTER LINE 

A I 3 IN (3 3 cm) FROM TOP 

▲ “ BOTTOM 



HORIZONTAL DISTANCE ( x) / WIDTH ( X) 


Figure 38. Exit plane velocity distribution along horizontal lines 
for dump diffuser; AR = k. 
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6-3 Cusp Diffuser 

The Inlet and diffuser lengths for the cusp diffuser were chosen 
to be equal to the corresponding lengths for both the Griffith and the 
dump diffusers. The area ratios were 2.7 and 3-7, as discussed in section A. 3- 

No standing vortices were observed using the available suction 
capacity and therefore attached flow was not possible to achieve. The 
flow emerged from the inlet section in the form of a jet and no pressure 
recovery was obtained whatsoever. The flow pattern is shown in Figure 
39- In all test runs the pressure in the inlet section was higher than 
the pressure in the exit plane (atmospheric). Referring again to Figure 
39 ♦ for an inlet velocity of li. = 100 ft/sec (30.5 m/s), typical values 
of the static pressure at several locations were as follows: 

P q - P. = 2.6 inch (6.6 cm H^O) 

P. - P c - 1.2 inch H 2 0 (3.0 cm H 2 0) 

P - P. = 1.0 inch H_0 (2.5^ cm H o 0) 

c t 2 i 

P* - P „ =0.0 inch H o 0 (0.0 cm H_0) 

t atm L 4 

P is the stream pressure at the cusp location, P is the pres- 
C t 

sure at the vena contracta, P. is the pressure at the inlet section and 
P Q is the stagnation pressure. These figures show that no pressure 
recovery was achieved and the channel did not act as a diffuser. It is 
believed that a longer diffuser length is needed to achieve unseparated 
flow with two standing vortices than the length which was used (6.2 in 
(15.7 cm) ) . 

An exit velocity map is shown in Figure which shows regions of 
backflow near the diffuser walls. 
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SECTION VI I 


CONCLUSIONS 

7.1 Gri ff i th Di ffuser 

In a metastable operating condition, slot suction rates around 8% 
of the inlet through flow rate was required to prevent flow separation 
for AR = *t, provided that adequate s i de-wal 1 suction was applied. For 
AR = 3, the values of slot suction rates were about 25% lower than those 
required for the higher area ratio. For stable operation, slot suction 
rates of approximately 30% were required with the diffuser of AR = k, and 
28% with the diffuser of AR = 3. For nearly all the unseparated test 
runs, the diffuser effectiveness n was 100%. 

7.2 Dump Diffuser 

The location of the suction slot proved to be critical as regards 
the diffuser effectiveness n. Much higher values of diffuser effective- 
ness were observed with the slot located at the inlet corner in contrast 
to the exit corner. However, flow separation occurred in all test runs. 

7. 3 Cusp Di ffuser 

Using the same inlet and diffuser lengths as for the Griffith 
and dump diffusers, 1 in (2.5*0 cm for inlet and 6.2 in (15-7) cm for 
diffuser lengths, it was not possible to achieve unseparated flow 
with two standing vortices in the cusp diffuser. The flow behaved as a 
jet emerging from the inlet section and no pressure recovery was 
obtained. It is believed that unseparated flow with two standing vortices 
could be achieved using a diffuser length longer than the one used in 
this investigation. 
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APPENDIX A 


Outline of Surface - Source - Distribution Method 
In Solving Channel Flow Problems Using Cascade Theory 

In determining the flow field around non-lifting bodies, the surface - 
source - distribution method is used. The basic idea is to cover the 
body surface with a continuous distribution of sources (or sinks) and 
to adjust the strength of this distribution so that the total velocity 
normal to the body surface is either zero (impermeable wall) or a 
prescribed value (permeable wall). Once the distribution of sources is 
known, the flow field is completely determined. The aim now is to outline 
the method by which this distribution could be obtained. 

Consider the flow around a two-dimensional body as shown in Figure 41 . 
is the velocity vector which represents the stream velocity upstream 
of the body. Let U n represent the prescribed velocity normal to the 
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surface, for a solid surface U = 0. If 4* is the potential due to the 

n 

source distribution at a particular point p on the surface, then the 
velocity normal to the surface at point p must be made up of two parts 
according to the equation 


U = U • n + 

n 00 3n 


32 


where n is the outward normal vector 

In this equation * n represents the component of normal to the 

A 

surface at point p. The term — represents the contribution of the 

on 

potential due to the yet unknown source distribution, henceforth, called 
disturbance potential. If the source - distribution per unit area at any 
point q is denoted by o(q), then the potential at p due to a source o(q)oS 
at q is given by the equation: 


<5<f>(p) = o(q) j&S 

The distrubance potential at p due to the total source distribution over 
the entire body surface S is given by: 


<f>(p) 


‘ // a(q) r 


dS 


33 


The point p need not be on the surface, as shown in Figure 4l , since 
equation 33 is equally applicable if p is off the body surface. The 
fluid velocities in the three directions of a Cartesian coordinate 
system x, y, and z, are 34>(p)/3x, 3<|>(p)/3z respectively. The velocity 
vector at point p is given by: 

u = 2 ,r.d*(p ) „ 

-r- ** 

where i, j, and k are the unit vectors along the x, y, and z coordinates 
respecti vely. 
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35 


From equations 33 and 
grad <f>(p) = 



, equation 35 can be derived: 
a(q) grad (-£-) dS 


Using equation 35 the boundary condition expressed in equation 32 can be 
written as: 

U - (U -*n) = Jj a[ q) (-J-) dS 36 

n 00 2 dn r 

For a known body contour and a known free stream velocity U , equation 
36 can be solved using numerical integration to yield the source 
distribution cr(q). Once the source distribution a(q) is known, the 
flow field is determined. 

In solving equation 36 numerically, the body contour is 
approximated by line segments and assuming that the source density 
a(q) is constant along each segment, see Figure 4l . Letting the number 
of segments be N, and using summation instead of integration, equation 
36 can now be written as: 


- N 

U - (U • n) = E A o (q) 
n - q= , pq 


37 


3 1 

where A = — (— ) AS and AS is the area of the element considered, 
pq dn r 

Equation 36 represents one equation in N unknown source strengths, a(q), 
where q varies from 1 to N. To obtain N equations, the point p is 
allowed to change position from element 1 to N. Hence, we obtain a set 
of N simultaneous equations in N unknowns. An expression analogous to 
that of 36 can be written for the tangential components due to the 


disturbance potential, 


9<{>(p) 

at 


. Thus , 


94>(p) 

at 


N 

E 

q=l 


pq 


o(q) 
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In equations 37 and 38. A and B are respectively, the normal and 

pq pq 

tangential velocities, at point p due to a unit source strength at point 
q. The total velocity due to the disturbance potential is the resultant 
of the normal and tangential velocities, and is given by equation 39 


N 

grad 4>(p) = E ( A p q 5 + B pq j) 0 (p) 39 

In order to evaluate the coefficients A and B at point p due 

pq pq 

to the influence of a source distribution of unit strength along the 
segment at q, the concept of complex velocity is used for convenience. 
For two-dimensional flow, the complex velocity at point p, w pq » < s given 
by: 


W = U - i U 
pq x y 

where U and U are the velocity components along the x and y axes 
x y 

pectively. 


40 

res- 
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Figure 42 shows an enlarged sketch of the segment at point p,< the rest of 
the body surface is shown dotted. The normal and tangential velocities 
at p due to a unit source distribution along the segment at q are shown; 

A and B respectively. The segment at p is shown to be inclined at 

pq pq 

an angle of X to the x-axis. From Figure 42 it can be shown that: 

U =B cos X - A sin X 

x pq pq 

and 

U =B sinX+A cos X 

y pq pq 


It is obvious from equations 40 and 4l that once the complex velocity 



is known, the coefficients A and B can be determined usinq equation 40 

pq pq 

To evaluate W , consider the complex veloci-ty at point p due to 
pq 

a source distribution of unit strength per unit length along the segment 
q. Figure ^3 shows an enlarged sketch of the segment q with point p 



/ 

/ 


Figure 43. Segment q with point p outside body surface 
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drawn outside the body surface. Let be the complex coordinates of 

point p: Z = x + iY ; and c be the complex coord i nates of point q; 

K K P P P q 

t = £ + in . It can be shown that the compl ex potent i a 1 F (defined as 

q q q 

<j> + iip) at point p due to a unit source located at ?(s) is given by: 

F = fe ln [Z p ‘ ds)] *»2- 

The complex velocity at point p is given by 

f-irsz 'V *s 

The complex velocity at point p due to the segment q having a unit source 
strength distribution per unit length. This can be derived by integrating 
equation b 3 along the length of the segment q from to ^ > see 
Figure 


pq 2 tt 


S * 


segment q 


In U p - C (s>] ds 


where ds is an infinitesimal length along the segment. The same Figure 
also illustrates the following relationship: 

ds e = dc, bS 

Replacing ds from equation b 5 into equation bb gives: 


2q 3T- (z p ' 5) dt 

S lq P 


- i Aq 

~2tt ln ^ Z p " ^ 


, <Z p - hq> 

~2T~ (f - cT 

P 2q 


The flow channel of the two-dimensional diffuser may be analyzed as 
a cascade problem, see Figure bb . The flow and the body - surface - source 
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density is identical for each body. Equation 46 is no longer valid 

for the cascade problem. It can be shown that the new expression for 

W is given by equation 47. 

pq 

t 



Figure 44. Flow channel represented by a cascade of bodies 


sink [ (tt/SP) U - e )] 
p iq 

sink LU/SP) (Z - s J] 
P 2q 


■ i Aq 


W 


pq 
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The numer 
solve N s 
in obtain 


cal computations were performed 
multaneous equations similar to 
ing the source distribution a(q) 


using a digital computer to 
equation 36 which result 
and hence the flow field. 
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APPENDIX B 


Derivation Of The Expression For The Diffuser Effectiveness 


The basic definition of the diffuser effectiveness is defined as 
the ratio between the actual and the ideal pressure recovery. 

P - P. 

n = (P - P.) 7 

e i ideal 

The ideal pressure recovery is derived from Bernoulli's equation as 


f o 1 1 ows : 


(P - P.) . . = | (u 2 - u 2 ) 

e i ideal 2 i e 


With no suction, U will be: 


Wi th suction, U wi 1 1 be: 
e 


U e=“AR U i 


% suction] 


Substituting this value of U e into equation 49 we get: 


m - % suction V 

(P - P.) . , , = (P - P.) 1 - 5 1 

e i ideal o i \ AR / 

* s 

1 2 

where (P q - P.) was substituted for y p IL in equation 51. Substituting 
this value of (P^ - Pj)j c j ea | ' nto equation 48, the final expression for the 
diffuser effectiveness is obtained 


P - P. 
e i 


(P_ " 


% suction 
100 
AR 


64 



FINAL REPORT DISTRIBUTION LIST FOR NASA CR-121024, Grant 41-001-031 


1. NASA - Lewis Research Center 

21000 Brookpark Road 
Cleveland, Ohio 44135 
Attention: Report Control Office 

Technology Utilization 
L i brary 

Fluid Systems Components Division 
W.L. Stewart 
L. Schopen 
J.B. Esgar 
W.T. Olson 
R.A. Rudey 
J.F. Dugan, Jr. 

Seymour Lieblein 
R.E. Jones 
Jack Grobman 
Lt . Col . G . S. Weden 
James A. Albers 
A.J. J.uhasz 

2. Mr. L. Magitti 
Naval Air Propulsion Test Center 
Bldg. 600 

Philadelphia Naval Base 
Philadelphia, Pennsylvania 19112 1 

3. Aerospace Research Laboratory 

Wr ight-Patterson AFB, Ohio 45433 

Attention: Dr. R.G. Dunn 1 

4. NASA Scientific and Technical Information Facility 
P.0. Box 33 

College Park, Maryland 20740 
Attention: NASA Representative 

RQT-2448 6 

5. FAA Headquarters 

800 Independence Avenue, S.W. 

Washington, D.C. 20533 

Attention: William Westfield 1 

Library 1 


MS 5-5 1 

3-19 1 

60-3 2 

5-3 1 

77-2 1 

500-206 1 

60-4 1 

3-16 1 

60-6 1 

86-1 1 

501-5 1 

60-6 1 

60-6 1 

500- 317 1 

501 - 2 1 

60-4 10 


65 



6. NASA Headquarters 

600 Independence Avenue, S.W. 

Washington, D.C. 20546 
Attention: N.F. Rekos (RAP) 

W.H. Roudebush (RAA) 

7. Department of the Army 

U.S. Army Aviation Material Laboratory 
Propulsion Division (SAUFE-PP) 

Fort Eustis, Virginia 2360*1 
Attention: J. White 

E.T. Johnson 

8. United Aircraft of Canada, Ltd. 

P . 0 . Box 

Lonquenil, Quebec, Canada 
Attention: Miss Mary Cullen 

9. Air Force Office of Scientific Research 
1 400 Wilson Boulevard 

Arlington, Virginia 22209 
Attention: SREP 

10. Defense Documentation Center (DDC) 

Cameron Station 

5010 Duke Street 
Alexandria, Virginia 22314 

11. Department of the Navy 

Bureau of Naval Weapons 
Washington, D.C. 20025 
Attention: Robert Brown, RAPP14 

12. Department of the Navy 

Bureau of Ships 
Washington, D.C. 20360 
Attention: G.L. Graves 

13. NASA-Langley Research Center 
Langley Station 

Technical Library 
Hampton, Virginia 23365 
Attention: Mark R. Nichols 

John V. Becker 

Richard J. Margajon MS 404 


1 

1 


1 

1 


1 


1 


1 


1 


1 


I 

1 

1 


66 



14. United States Air Force 
Aero Propulsion Laboratory 
Area B, Bldg. 18D 

Wr ight-Patterson A.F.B. 

Dayton, Ohio 45*03 
Attention: Robert E. Henderson 

15. United Aircraft Corporation 
Pratt S Whitney Aircraft Division 
400 Main Street 

East Hartford, Connecticut 06108 
Attention: G. Andreini 

L i brary 
R. Marshall 

16. United Aircraft Research 

East Hartford, Connecticut 
Attention: Library 

17. Detroit Diesel Allison Division 
Department 8894, Plant 8 

P.0. Box 894 

Indianapolis, Indiana 46206 
Attention: J.N. Barney 

G.E. Holbrook 
Library 

18. Northern Research & Engineering Corp. 
219 Vassar Street 

Cambridge, Massachusetts 02139 
Attention: K. Ginwala 

19- General Electric Company 
Flight Propulsion Division 


C i nc i nnat i , 

Ohio 45215 


Attention : 

J.S. McBridge 

H-44 


F. Burggraf 

H-32 


C. Danforth 

H-32 


Technical Information 

Center N-32 


D. Bahr 

20. General Electric Company 
1000 Western Avenue 
West Lynn, Massachusetts 01905 
Attention: Dr. C.W. Smith 

Library Bui lding 


l 


1 

1 

1 


1 


1 

1 

1 


1 


1 

1 

1 

1 

1 

1 


2-40M 


1 



21. Curt i ss-Wr ight Corporation 

Wright Aeronautical Division 
Wood-Ridge, New Jersey 07075 
Attention: D. Wagner 

W. Walker 

22. Air Research Manufacturing Company 
k02 South 36th Street 

Phoenix, Arizona 85034 
Attention: Robert 0. Bullock 

23. Air Research Manufacturing Company 
9851 Sepulveda Boulevard 

Los Angeles, California 90009 
Attention: Mr. L.C. Wright 

2k. AVC0 Corporation 
Lycoming Division 
550 South Main Street 
Stratford, Connecticut 
Attention: Claus W. Bolton 

Charles Kuintzle 

25. Continental Aviation & Engineering Corporation 
12700 Kercheval 

Detroit, Michigan 48215 
Attention: Eli H. Bernstein 

Howard C. Walch 

26. International Harvester Company 
Solar Divison 

2200 Pacific Highway 

San Diego, California 92112 

Attention: P.A. Pitt 

27. Goodyear Atomic Corporation 
Box 628 

Piketon, Ohio 

Attention: C.0. Langebrake 

28. George Derderian AIR 53622 B 
Department of the Navy 
Bureau of the Navy 
Washington, D.C. 20360 


1 

1 


1 


1 


1 

1 


1 

1 


1 


1 


1 


68 



29. The Boeing Company 
Commercial Airplane Division 
P.0. Box 3991 

Seattle, Washington 98124 

Attention: G.J. Schott MS 80-66 l 

30. The Boeing Company 

Missile and Information Systems Division 


224 N. Wilkinson Street 
Dayton, Ohio 45402 

Attention: Warren K. Thorson 1 

31. Aerojet-General Corporation 
Sacramento, California 95809 

Attention: M.S. Nylin 1 

Library 1 

32. Cornell Aeronautical Laboratory 
4455 Genessee Street 

Buffalo, New York 14221 1 

33* Marquardt Corporation 
16555 Saticoy Street 

Van Nuys , Cali forn i a 1 

34. Thompson Ramo Wooldridge 
23555 Euclid Avenue 

Cleveland, Ohio 1 

35- Aro, Incorporated 

Arnold Air Force Station 

Tennessee 1 


36 . The Johns Hopkins University 
Applied Physics Laboratory 
Silver Spring, Maryland 20910 
8621 Georgia Ave. 

Attn: W.B. Shippen 1 

Dr. G. Dugger 1 

37. Cummings Engine Co. 

Cummings Technical Center 
1900 McKinley Avenue 
Columbus, Indiana 47201 
Attention: Curt Dasbach 

Mail Code 50142 1 


69 


38. Garrett/Ai Research Co. 

402 South 36th Street 
Phoenix, Arizona 8503** 

Attention: John M. Haasis 1 

39 . Pratt £ Whitney Aircraft 

Florida Research £ Development Center 
Box 2691 

West Palm Beach, Florida 33**02 


Attention: J. Chamberlain 1 

J. Dykslag 1 

J. Shadowen 1 

G. Lewis 1 


40. Professor A.H. Lefebre 

The Cranfield Institute of Technology 
Cranfield, Bedford 

Great Britain 1 

41. Aerojet General Corporation 
Sacramento Facility 

P.0. Box 15847 

Sacramento, California 95813 


Attention: C.E. Tedmon 1 

Dave Kors 1 

42. The University of Toledo 
Toledo, Ohio 43606 

Attention: Dr. Duen-Ten Jeng 1 

Dr. Kenneth DeWitt 1 

43. Eaton Yale and Towne Research Center 
26201 Northwestern Highway 

Southfield, Michigan 48075 1 


44. Rocketdyne 

North American Rockwell 
6630 Canoga Avenue 
Canoga Park, California 91304 
Attention: S.D. Clapp 

Manager 

Propulsion Technology Research Division 1 

45. Dept, of Mechanical Engineering 
201 Engineering Building 
Michigan State University 

E. Lansing, Michigan 48823 

Attention: Dr. Merle C. Potter 2 


70 



46. Department of Engineering 
Leicester University 
Leicester, Great Britain 

Attention: Dr. David J. Cockrell 1 

47. Von Karman Institute for Fluid Dynamics 
Chaussee de Waterloo 72 

1 640 Rhode-St-Gense (Beligique) 

Attention: Professor J. Chaurvin 1 

48. Simon Laboratories 
University of Manchester 
Manchester, Great Britain 

Attention: Dr. John Turner 1 

49. Engineering Department 
Salford University 
Manchester, Great Britain 

Attention: Dr. Wilcox 1 

50. Engineering Department 

University of Liverpool 
Liverpool, Great Britain 
Attention: Professor J.C. Gibbings 

51. Gas Turbine Department 
General Electric Co. 

Schenectady, New York 

Attention: R.L. Hendrickson 1 

52. School of Mechanical and Aerospace Engineering 
Upson Hall 

Cornel 1 Un i vers i ty 
Ithaca, New York 

Attention: Professor E.R. Resler 1 

Professor D.G. Shaphard 1 

53. School of Engineering 
University of Auckland 
Privatebag, Auckland 
New Zealand 

Attention: Dr. G.L. Bowen . 1 


71 



54. Technical Center 
Cummings Engineering Co. 

Columbus, Indiana 

Attention: Dr. A.C. Chu 1 

55. Mechanical Engineering Department 
Iowa State University 

Ames, Iowa 50010 

Attention: Professor George K. Serovy 1 

56. Mechanical and Aerospace Engineering Dept. 

University of Tennessee 

Knoxville, Tennessee 

Attention: Professor Mancil Milligan 1 

57- Department of Mechanical Engineering 

Virginia Polytechnic Inst. & State Univ. 

Blacksburg, Virginia 24061 

Attention: Professor Felix J. Pierce 1 

58. Department of Mechanical Engineering 
Purdue Un i vers i ty 

Lafayette, Indiana 

Attention: Dr. A.T. McDonald 1 

59. McDona 1 d-Doug las Aircraft Co. 

St. Louis, Missouri 63 1 66 

Attention: Dr. E.D. Spong 1 

Dept. 243 

60. Aerospace Group System Division 
Hughes Aircraft Co. 

Canoga Park, California 91304 

Attention: Mr. J.L. Raymond 1 

61. Lehigh University 
Bethlehan, Pennsylvania 18015 

Attention: Dr. Donald 0. Rockwell 1 


72 



